Optical Observations of the Binary Pulsar System PSR B1718— 19: 
Implications for Tidal Circularization 



M. H. van Kerkwyk^-^, V. M. Kaspi^ A. R. Klemola^ S. R. Kulkarni^ A. G. Lyne^ 

and 

"qv^ ' D. Van Buren'^ 



< 



in 



ABSTRACT 



We report on Keck and Hubble Space Telescope optical observations of the eclipsing 
i binary pulsar system PSR B1718— 19, in the direction of the globular cluster NGC 6342. 

These reveal a faint star (mF702W = 25.21 ± 0.07; Vega system) within the pulsar's O'.'S 
radius positional error circle. This may be the companion. If it is a main-sequence 
star in the cluster, it has radius Rq — 0.3 Rq, temperature Teg ~ 3600 K, and mass 
Mq ~ 0.3 Mq. In many formation models, however, the pulsar (spun up by accretion 
00 ' or newly formed) and its companion are initially in an eccentric orbit. If so, for tidal 

^ ■ circularization to have produced the present-day highly circular orbit, a large stellar 

'. radius is required, i.e., the star must be bloated. Using constraints on the radius and 

"Th I temperature from the Roche and Hayashi limits, we infer from our observations that 

Ph' Rc ^ 0.44 /?0 and Tcs ^ 3300 K. Even for the largest radii, the required efficiency of 

Q ■ tidal dissipation is larger than expected for some prescriptions. 

■ Subject headings: binaries: close — pulsars: individual (PSR B1718— 19) — stars: evo- 

lution 
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1. Introduction 

PSR B1718— 19 is a 1-s radio pul sar in the directio n 
of the globular cluster NGC 6342 ( |Lyne et al. 1993|) . 
It is in a 6.2 h circular binary orbit with a compan- 
ion that has mass Mq > 0.11 Mq (assuming a pulsar 
mass AfpsR — 1.35 Mq). At radio frequencies around 
400-600 MHz, irregular eclipses of the pulsar occur 
around superior conjunction, indicating the presence 
of material around the companion. The eclipsing ma- 
terial must be tenuous since the eclipses are absent 
at higher frequencies. The properties of the pulsar 
itself are not special; from the spin period and its 
derivative, one infers a dipole magnetic field strength 
B ~ 3.2xl0i9(PP)i/2 ~ 1.5 X 10^2 G and a spin-down 
ageTsd = -P/2P = 10Myr. 

The system appears peculiar for several reasons 
( Lyne et al. 1993[ ): (i) supernova activity in NGC 6342 
should have ceased long ago, hence one does not 
expect it to contain an apparently recently formed 
pulsar; (ii) all other known binary pulsars in globu- 
lar clusters have millisecond periods and low mag- 
netic fields, evidence that they were "recycled" or 
spun up by accretion from a binary companion, the 
field presumably having decayed in the process (for 
a review, see Van den Heuvel 1995); (iii) unlike in 



other eclipsing pulsars, the rotational energy loss of 
PSR B1718— 19 is many orders of magnitude below 
the critical flux necessary to drive a wind from the 
companion, so the origin of the eclipsing material is 
unclear. 
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have focused on the three peculiarities discussed 
above, mostly in the context of the two favored for- 
mation scenarios: forming a new pulsar by accretion- 
induced collapse of a white dwarf, and recycling an 
old pulsar in a close encounter with other stars in the 
core of the globular cluster. We believe, however, an- 
other important issue is that (iv) the current orbit is 
near-circular (e <^ 0.005), while the formation should 
have left the pulsar and its companion in an eccentric 
orbit. An initially eccentric orbit is not only expected 
in all formation models, but also indicated by the fact 
that the system currently is not in the core of the clus- 
ter, as would be expected given its mass, but offset by 
2.'4; apparently, it received a kick, which should have 
made the orbit eccentric as well. 

It has been noted by Verbunt ( |1994[ ) that if circu- 



larization occurred, the energy dissipated in the com- 
panion would have been of order the binding energy of 
the star. Thus, the star might have become bloated or 
even have been partially destroyed. Little attention, 
however, has been given to the question of whether 
the tidal circularization efhciency is sufhcient for cir- 
cularization to have happened within the spin-down 
age. 

The above puzzles motivated us to try to identify 
the counterpart of PSR B1718— 19. In this article, 
we report on Jodrell Bank, Very Large Array (VLA), 
Keck and Hubble Space Telescope (HST) observations 
of PSR B1718-19. The improved position of the pul- 
sar from the radio observations and the results of the 
optical imaging observations are reported and sum- 
marized in §1]. Given the crowded field, establishing 
a precise tie between the radio and optical observa- 
tions is crucial. A detailed account of our astrometry 
is given in The photometry of a candidate object 
is also reported in this section. The constraints on the 
basic parameters of the companion derived from op- 
tical and pulsar timing observations are discussed in 
In §1^, we discuss the implications of our results, 
focusing in particular on the circularization. 

2. Observations 

The position of the pulsar was first determined 
from radio timing measurement s using the 76-m 
Lovell telescope at Jodrell Bank ( Lyne et al. 1993| ). 
We derived an improved timing position from data 
obtained between 1994.6 and 1998.8: aj2ooo = 
17h21'°01!53 and 5j2ooo = -19°36'36", with uncer- 
tainties of 0!04 and 6", respectively. This position has 
a relatively large uncertainty in declination because of 
the source's proximity to the plane of the ecliptic. An 
image of the field was obtained on 11 November 1992 
using the VLA in the A-configuration. It reveals a 
source within the timing error box with a flux den- 
sity of 0.3 mjy, consistent with the flux of the pulsar 
at the observing frequency of around 1400 MHz. The 
position of this source is aj2ooo — 17^21™01!54 and 
(^.Kooo = — 19°36'36'.'6, with an micertainty of 0'.'2 in 
each coordinate. 

The PSR B1718-19 field was observed on 4 June 
1995 using the Low-Resolution Imaging Spectrome- 
ter (LRIS; |Oke et al. 1995| ) on the 10-m Keck tele- 
scope. Two series of exposures with integration times 
ranging from 10 s to 240 s were taken through several 
filters. The seeing was ~ 0'.'9 during the first series. 
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Fig. 1. — Optical images of the field of 
PSR B1718-19. The top panel shows part of the 
LRIS I-band image. The exposure time was 240 s 
and the seeing 0'.'9. The limit to the flux of any 
star within the (X'S radius (95% confidence) error 
circle centered on the VLA position is / > 24.0 
(95% confidence). The bottom panel shows part of 
the Planetary Camera image, taken through the the 
F702W filter. The image shown is the median of all 
twelve 700 s exposures. The candidate counterpart 
has TOF702W = 25.21 ± 0.07. 



and ~ 1'.'2 during the second. Our best-seeing I-band 
0.8/xm) image is shown in Fig. ^ These images 
showed that there was no object in the VLA error 
box, and that with 95% confidence / > 24.0 mag. 
Furthermore, they showed that a brighter star was 
so close that it would be hard to make any further 
progress from the ground. 



The Wide Field Planetary Camera 2 (WFPC2) on 
board the HST was used to observe the field through 
the F702W filter in four contiguous HST orbits, from 
14:39 to 20:07 UT on 8 March 1997. In each spacecraft 
orbit, three 700-s exposures were taken, offset in both 
X and Y by 0, 3, and 6 PC pixels in order to miti- 
gate the effects of hot pixels and imperfect flat fields. 
The pointing was such that the field around the pul- 
sar is on a clean spot on the CCD of the Planetary 
Camera (PC). 

The reduction and analysis started from the 
pipeline-calibrated PC images (Holtzman et al 



1995a). We first obtained median images for each 
of the three observing positions. The median image 
for the middle offset was used for the astrometry de- 
scribed below. Next, for each of the three sets, we 
found pixels hit by cosmic rays in the individual im- 
ages by comparison with the median image; pixels 
with values more than 5cr above the median were re- 
placed with the median (here, a is an estimate of the 
expected uncertainty based on the median value). Fi- 
nally, the images were registered by applying integer 
pixel shifts and co-added to form a grand average. 
This average is free of cosmic rays, but not of hot 
pixels and other chip defects. We used it for the pho- 
tometry nonetheless, as our candidate (see below) and 
the other stars we selected were on clean parts of the 
chip. For display purposes, however, we have used 
the median of all registered images in Fig. 

3. Astrometry and Photometry 

The astrometry of the PSR B1718-19 field was 
carried out in four stages. First, 23 stars from the 



ACT Reference Catalog ( [Urban, Corbin, fc Wycofl 
1998b) were used to derive an astrometric solution 
for a plate taken at epoch 1988.37 on Kodak 103a-G 
emulsion using the yellow lens of the 0.5-m Carnegie 
double astrograph at Lick Observatory. The model 
employed five terms in each coordinate (proportional 
to 1, X, y, xy, and {x'^,y'^) in {x,y)). The inferred 
rms error for a single star is 0'.'16 in each coordinate, 
and the zero-point uncertainty in the solution 0'.'04. 
We note that the errors are larger by about 50% than 
expected based on the measurement and ACT coor- 
dinate uncertainties. Probably, this is because some 
stars are blended with fainter objects; thus, it should 
induce no systematic error. 

Second, the plate solution was used to calculate 
right ascension and declination for 30 fainter, rela- 
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tively isolated stars in common with the short LRIS 
R-band exposure. For these stars, positions were also 
measured on the LRIS frame, and corrected for in- 
strumental distortion using a bi-cubic transformation 
determined by J. Cohen (1997, private communica- 
tion). Solving for offset, scale and rotation, the in- 
ferred rms single-star error is 0'.'26 in each coordinate, 
consistent with the expected measurement errors on 
the astrograph plate for these relatively faint stars. 
The zero-point uncertainty in the solution is 0'.'05. 

Third, using 23 fainter stars on a smaller part 
of the LRIS image around the VLA position of 
PSR B1718— 19, the solution was transferred to one 
of the 240 s LRIS R-band images, solving only for the 
offset between the two exposures. The inferred rms 
error is O'.'Oll in each coordinate, and the zero-point 
uncertainty in the solution 0'.'002. 

Finally, 37 stars were used to tie the astrometry to 
the WFPC2 PC image (14 of these were used in the 
previous stage as well). The WFPC2 positions were 
corrected for instrumental distortion using the cubic 
transformation given by Holtzman et al. ( 1995a ). Off- 
set, scale and rotation were left free in the solution^. 
The inferred rms error was 0'.'017 in each coordinate, 
and the zero-point uncertainty in the solution 0'.'003. 

From the numbers above, the total uncertainty in 
the, zero points of positions due to measurement er- 



ror « is 0"06 in each coordinate , We expect S3rstcmatic 



effects to be smaller than this; e.g., no large system- 
atic offset as a function of magnitude is seen on yel- 
low astrograph plates used for other applications. A 
possible additional uncertainty is the extent to which 
our ACT-based solution is on the same system as the 
VLA position of PSR B1718-19, i.e., the Interna- 
tional Celestial Refer ence Syste m (ICRS). Th e ACT 
co mbines th e Tycho ( |ESA 1997D and AC2000 ( [Urban 
et ;d. 1998a) catalogs to determine positions and ac- 
curate proper motions. The Tycho positions are on 
the ICRS to within 0.6 mas ( |ESA 1997| ), and any off- 
set from the ICRS at the plate epoch, which is close 
to the Tycho epoch (1991.25), should be small as well. 

In summary, we expect the uncertainties in the 
astrometry to be dominated by measurement er- 



* The fitted scale is close to the value given by Holtzman et 
al. ( |l995a| ), as is the ofli'set between the fitted position angle 
and the position angle of the HST V3 axis listed in the image 
header. If we fix the values, the inferred position of our pro- 
posed counterpart changes by 0''006 in RA and — O'.'OOS in Dec. 
Thus, a perhaps more conservative estimate of the zero-point 
uncertainty is about 0''01 in each coordinate. 



rors, 0'.'06 in each coordinate for the optical posi- 
tion, and 0'.'2 for the radio position. Combining 
the two in quadrature, the 95% confidence error ra- 
dius is [-21og(l - 0.95)(0'.'22 -|- 0'.'062)]i/2 ^ o'.'5. 
Within this radius, there is one faint object (Fig. 
at aj2ooo = 17h21™01!549 ± 0';004 and <5j2ooo = 
-19°36'36'.'76 ± 0'.'06. This may be the optical coun- 
terpart of PSR B1718-19. 

We measured magnitudes for the candidate coun- 
terpart of PSR B1718— 19 and for a number of other 
stars in the field following the prescriptions of Holtz- 
man et al. ( |1995h| ) and Baggett et al. ( |1997D . We per- 
formed aperture photometry for a range of different 
radii, and used some thirty brighter stars in the frame 
to determine aperture corrections relative to the stan- 
dard 0'.'5 (llpix) radius aperture. For the candidate, 
the best signal-to-noise ratio is for relatively small 
apertures, with radii between 1.5 and 3pix. From 
these, we infer a count rate for the 0'.'5 radius aper- 
ture of 0.070 ± 0.005 DNs-i (we find consistent val- 
ues for the larger apertures; here, 1 DN corresponds 
to about 7 electrons, given the gain used for our ob- 
servations). This count rate corresponds to a mag- 
nitude mF702w = 25.21 ± 0.07 in the Vega system 
(using TOF702W = 22.428 for a count rate of 1 DNs~^ 
in the PC, and applying a 0.10 mag aperture correc- 
tion from 0'.'5 radius to "nominal infinity" ; Baggett et 



al. 1997). Magnitudes for the individual frames show 



a standard deviation of 0.25 mag around the average, 
similar to what is found for other stars at this bright- 
ness level. We do not find modulation at the 6^2 or- 
bital period of PSR B1718— 19 (which is well covered 
by our observations), although the limit on the mod- 
ulation amplitude is not very restrictive: < 0.3 mag 
at 95% confidence. 

At or above the brightness level of the candidate 
counterpart, there is about one object per four square 
arcseconds in this field. Thus, there is a probability of 
about one in five of finding an object within the 95% 
confidence error radius by chance. If it is a chance 
coincidence, the real counterpart must be substan- 
tially fainter. We derive a 95% confidence limit of 
"T'F702W = 27.2 for any other object in the error cir- 
cle (using the observed noise in the sky near the VLA 
position, of Csky = 5.6 x 10~'*DNpix^^ s^^, and the 
fact that within the 0'.'5 radius error circle there are 
about 400 resolution elements). 
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Observational Constraints 



the binary orbit (Lyne et al. 1993). From tlic mass 
function, assuming MpgR — 1.35 Mq (Thorsett & 



tern by the observations. We first summarize the 
constraints on the age of the system and the mass 
of the companion set by radio observations, and then 
discuss the constraints on the companion radius and 
temperature from our T/iST detection and Keck hmit. 
We will assume that we detected the companion, and 
that the system is located in NGC 6342. 

4.1. System Age 

Clearly, something happened to the system re- 
cently (as compared to the cluster age). One indica- 
tion is the short inferred spin-down age of the pulsar, 
Tsd = P/2P — lOMyr. This is the time required for a 
dipole rotating in vacuo to spin down to the present- 
day spin period from an infinitely fast rotation rate. 
It is thus an upper limit to the true age, unless other 
mechanisms influenced the spin period (e.g., transient 
spin-up by accretion or a braking index very different 
from that predicted by dipole emission). 

Another indication of a recent event is that the 
system is offset from the cluster core. Unperturbed, 
a system as massive as this should have settled in the 
core long ago due to mass segregation (in ^0.5 Gyr 
given its present position). Indeed, the progenitor(s) 
of this system, which must have been massive as well, 
should have resided in the core. It seems natural to 
argue that it was a single event that brought the sys- 
tem to its present state and kicked it out of the core. 

We note that it would require fine-tuning for the 
kick to have resulted in the system remaining in a 
cluster orbit as wide as is indicated by the 2.'4 off- 
set (the half-mass radius is 0.'9; for a discussion, see. 



yVe now evaluate the conscramrs set on the sys- Chakrabarty 1999| ), one infers Mc > 0.11 Mq; fur 



e.g., Phinney 1992). If instead the system is unbound. 



the kick needs to have happened <^ 1 Myr ago, which 
would imply that the event left the pulsar with a spin 
period only slightly shorter than the present one. This 
would be consistent with scenarios in which the pul- 
sar was spun up by accretion: for Eddington-limited 
accretion, the equ ilibrium spin period is close to 1 s 
( Lyne et al. 1993| ). We conclude that the system has 
been in its present state for 10 Myr at most. 

4.2. Companion Mass 

Radio pulse timing of PSR B1718— 19 has yielded 
the mass function, /(Mc, Mpsr) = sin^ V(^C + 
Mpsn)^ = 0.000706 Mq, where i is the inclination of 



thermore, there is a 95% a priori probability that 
i > 18° and hence that Mq < 0.43 Mq. 

The radio eclipses allow a direct constraint on i, 
albeit in a model-dependent way. Burderi & King 
(1994) and Thorsett (1995) calculated the expected 
attenuation of the radio flux at different frequen- 
cies for a simple constant-velocity, spherically sym- 
metric wind, and found that they could reproduce 
the observations of Lyne et al. (1993) for i > 30°. 
This would correspond to Mq < 0.25 Mq. Burderi 
& King argue that almost certainly i > 20°, im- 
plying Mc < 0.35 Mq. In summary, most likely 
0.11 ^Mc^ 0.35 Mq. 



4.3. Companion Radius and Temperature 

The interpretation of the apparent F702W mag- 
nitude depends on the companion radius Rc, effec- 
tive temperature T^g, and (to a lesser extent) surface 
gravity log 5, as well as the distance and reddening 
to NGC 6342. To determine the constraint set by 
our measurement, we used parameters for NGC 6342 
from the May 1997 edition of the catalog of globu- 
lar clusters (Harris 1996): reddening Eb-v — 0.44, 
distance modulus (m — M)v = 16.15, and metallic- 
ity relative to solar [Fe/H] = —0.65. The distance 



scale used in the catalog is similar to the Hipparcos- 
based one (e.g., that of Chaboyer et al. |1998| would 
give (m — M)v — 16.22). The reddening corre- 
sponds to Ay = 1.36, Ar — 1.02, A l — 0.65 (using 



the extinction curve of Mathis 1990 ), and one infers 
(m - M)o = 14.79 and d = 9.1 kpc. We also used the 
evoluti onary tracks for [Fe/H] = —0.5 stars of Baraffe 
et al. ( 1998| ), to relate temperatures and radii to ab- 
solute magnitudes in various standard bands. 

With these data in hand, we proceeded as follows. 
First, we looked up temperatures and corresponding 
radii for the main sequence (at age 1 Gyr, appropri- 
ate fo r a metal-rich globular cluster; Salaris fc Weiss 
199S). Second, we used My and Mn for these stars 
to calculate Vms and i?ins for the distance and red- 
dening of NGC 6342. Third, for each { V - R) ^^, we 



used the calibration of Holtzman et al. ( 1995b ; Eq. 9, 
Table 10) to infer the expected i?F702W magnitude 
corresponding to the observed F702W count rate. 
Fourth, we derived the radii required to match the 
observations from the difference i?F702W ~ ^^ms- To 
estimate the uncertainty, we assumed a total uncer- 
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tainty of 0.3 mag in the magnitude difference, which 
we regard as a 95% confidence estimate (it is domi- 
nated by the estimated uncertainties in distance and 
reddening). Strictly speaking, one should redden the 
F702W magnitude and then calculate i?F702W using 
{V — R)o. Furthermore, one should consider the influ- 
ence of log g on the star's colors, and take into account 
the slight difference in metallicity. None of these cor- 
rections, however, is important at the present level of 
accuracy. 

Another constraint is set by our I-band non- 
detection. This corresponds to an upper limit to 
the mF702W ~ I color, and thus a lower limit to the 
effective temperature. We find T^g > 3150 K for 
mF702W — I > 1.4 (here, we have decreased the limit 
to the magnitude difference by 0.2 mag to account for 
uncertainties in mF702Wi in the color transformation, 
in the effects of changes in \ogg, and in the redden- 
ing). 

The observational constraints can be summarized 
graphically in a diagram of Rc versus the effective 
temperature Toff of the companion (see Fig. ^). To 
make further progress we need to have some knowl- 
edge of the companion's nature. We consider two 
cases: the companion is an ordinary main sequence 
star; and the companion is a bloated star. 

4-. 3.1. A Main Sequence Companion 

In Fig. H, the locus of radius versus effective tem- 
perature for main sequence stars of different masses 
is indicated (from Baraffe et al. 199 j ). From the in- 
tersection with the region allowed by our observation, 
one infers that if the companion were a main sequence 
star, it would need to have 0.26 ^ Rq <, 0.33 i?Q 
and 3580 ^ T^ff ^ 3680 K. Its mass wouldl)e 0.26 ^ 
Mq <^ 0.34 A/0, consistent with the constraints in- 



ferred dynamically and from eclipse modeling (§4.2). 

4.3.2. A Bloated Companion 

The companion does not necessarily have to be on 
the main sequence. Indeed, bloating might be ex- 
pected: most formation mechanisms produce a bi- 
nary that has substantial eccentricity initially, and 
the amount of energy that needs to be dissipated in 
order to circularize the orbit is a substantial fraction 
of th e binding energy of a low-mass starQ ( |Verbunt 
1994 ). Bloating may also be required, as alluded to 
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Fig. 2. — Constraints on the radius and temperature 
of the companion of PSR B1718— 19. The continuous 
curve with dashed curves next to it indicates the con- 
straint set by our measured F702W magnitude. The 
bold vertical line with the hashed region to the right 
indicates the limit to the temperature set by the limit 
to the mFY02W—I color. The continuous curve labeled 
with masses (in solar units) indicates the relation ex- 
pected for main sequence stars with [Fe/H] = —0.5; 
the bold section of this curve represents the allowed 



range if the star is on the main sequence (§4.3.1). 
As argued in the text, the companion may have be- 
come bloated due to irradiation or tidal dissipation. 
The dotted lines indicate pre-main-sequence tracks 
for various masses. These are close to the Hayashi 
limit, and stars in hydrostatic equilibrium can only 
be on or to the left of these tracks. The long-dashed 
curve connects the maximum (Roche) radii the com- 
panion could have if it were on such a track. The 
shaded region shows the full range of allowed param- 
eters (§ [1.3.2| ). 



in § |T| and discussed in § ^ for tidal dissipation to have 
circularized the orbit in the short time since forma- 
tion. 

We are not aware of detailed calculations of bloat- 
ing due to tidal dissipation for almost completely con- 



^The star cannot be large because of having evolved off the 
main-sequence, as a star at or past the turn-off mass would 



have a luminosity much higher than the limit implied by our 
observations. Even if mass were lost, the thermal timescale for 
such a star likely is too long for the luminosity to have changed 
substantially. 
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vective, low-mass starsQ. Regardless of the expansion 
process, however, for a given mass an upper limit to 
the radius is set by the Roche radius, above which 
mass transfer would occur. Furthermore, for given 
mass and radius, a lower limit to the temperature is 



set by the Hayashi (1961) limit, below which a star 
cannot be in hydrostatic equilibrium. 

To delineate the constraint on the temperature, 
we can use pre-main-sequence tracks, since these fol- 
low the Hayashi limit closely. Indeed, it seems not 
unlikely that the companion is currently contract- 
ing along a similar track, since all processes that 
could have induced bloating (irradiation, tidal heat- 
ing) should have ceased to be operative (the pulsar's 
spin-down luminosity being very low and tidal dissi- 
pation having ceased as the orbit became circular). 
Of course, it is not clear that there has been enough 
time for the companion to adapt to a quasi-pre main 
sequence configuration. On the other hand, if the 
bloating is due to tidal heating, and if the energy was 
dissipated somewhere in the convective regions, con- 
vection would have distributed the energy throughout 
the star and expansion may well have been along a 
quasi-pre main sequence track as well. 

In Fig. H, dotted lines indicate the pre-main- 
sequence tracks of BarafFe et al. ( 1998| ), and the long- 
dashed line shows the Roche radii for the whole range 
of masses. Under our assumptions, the companion 
must lie between the latter line and the main se- 
quence, and at any given radius and temperature, 
an upper limit to its mass is set by the mass for 
which the pre-main-sequence track passes through 
that radius and temperatures. In addition, it must 
reproduce the observed F702W flux and be consistent 
with the mF702W — I limit. The allowed ranges are 
3300 ^ Teff < 3680 K and 0.44 ^ i?c > 0.26 i?0. The 
inferred masses are 0.11 < Mq < 0.34 M©, consistent 
with the constraints inferred dynamically and from 



eclipse modeling (§4.2) except for the very low-mass 
end {Mq ^ 0.12 A/q), at which total eclipses would 
be expected (for a 1.35 A/q neutron star), which are 
not observed. 

5. Discussion 

We now discuss the implications of our results. The 
assumption we continue to make is that the system 
is associated with NGC 6342 and that we have de- 



tected the counterpart. The constraints we regard as 
most important are: (i) the system was brought to 
its present state <^ 10 Myr ago; (ii) the current orbit 
is nearly circular (e < 0.005); and (iii) the companion 
has 0.26 <, Rc <, 0.44 i?© and 3680 ^ T^a ^ 3300 K. 
In these ranges, the left-hand limits correspond to 
the case where the companion is a ~ 0.3 Mq main- 
sequence star, while the right-hand limits correspond 
to the case where the companion is a ~0.13A'/q star 
which has been maximally bloated and is currently 
contracting along a pre-main-sequence track. 

Below, we first briefly review the formation models 
that have been suggested, in order to set the stage for 
a discussion of how the system could have been cir- 
cularized. We end by noting briefly what constraints 
one could set if the star we detected is not the coun- 
terpart of PSR B1718— 19, or if the system is not in 
NGC 6342. 

5.1. Formation Models 

Two models have been suggested to explain the 
origin of the PSR B 1718— 19 system. One is that the 
neutron star was formed early in the life of the globu- 



lar cluster ( 


Lyne et al. 1993; 


Wijcrs & Paczynski 1993 


Zwittcr 1993; Ergma et al. 


1996 


). It had stopped be- 



ing a radio pulsar and was dormant until it had a close 
interaction with a star or binary, ~ 10 Myr ago. Dur- 
ing that interaction, some mass was lost from a nor- 
mal star, part of which was accreted by the neutron 
star. As a consequence, it was spun up sufficiently for 
the radio mechanism to become active again, yet did 
not undergo a long phase of mass transfer and thus 
kept a large magnetic field. Similar close encounters 
have been invoked to explain the presence of two other 
long-period (but single) pulsars in globular clusters: 
PSR B1745-20 in NGC 6440 and PSR B1820-30B 



in NGC 6624 (Lyne, Manchester, fc D'Amico 1996) 



To give the system the velocity required by the 
present location far outside the cluster core and con- 
serve momentum, the putative interaction must have 
involved at least one other object. The system's ini- 
tial orbit is expected to have been highly eccentric. 



^"See Podsiadlowski ( |l996[ ) for calculations for somewhat more 
massive stars, for which only the outer layers are convective. 



e ^ 0.8 (see, e.g, [Phinncy 1992| ). 

The second model that has been considered for the 
origin of the PSR B1718— 19 system is that the neu- 
tron star was formed via accretion-induced collapse 
(AIC ) of a white dwarf (Fync et al. 1993| ; |Ergma| 
199^ ; IWijers fc Paczynski 1993| ; Ergma et al. |1996|) . 



The kick imparted to the neutron star during AIC, 
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whether intrinsic or due to mass loss, must have been 
relatively small, since otherwise the systemic veloc- 
ity would have been so large that the system would 
have left the cluster long ago. For a kick of order 



the 20 km s ^ escape velocity from the core ( Webbink 
19^, an initial eccentricity of ~0.2 is expected. 

Given such a kick, the increase in the orbital sepa- 
ration should have been small, about a factor 1.25. 
Thus, the companion, which must have filled its 
Roche lobe before AIC for mass transfer to occur, 
should still be close to filling its Roche lobe. This 
implies that it cannot have been a Roche-lobe filling 
main-sequence star, as this would require a mass of 
^0.7 Mq and hence a much brighter optical counter- 
part than we observe (see Fig. ft has been sug- 
gested that the star was bloated already before AIC 
due to irradiati on from the primary duri ng the mass 
transfer phase ( Ergma 1993 ; Ergma et al. 1996 ). ft is 
not clear whether strong bloating i s possible by irra - 
diating only one side of a star (e.g., King et al. 1996 ), 
but if it happened, the star should still be bloated, 
because the thermal time scale is much longer than 
the pulsar characteristic age. 

5.2. Circularization 

In both formation scenarios that have been sug- 
gested, the initial binary orbit is expected to be ec- 
centric. Indeed, in general any event that imparted 
a systemic velocity large enough for the system to 
(almost) escape the cluster will most likely have left 
the binary orbit eccentric. The current tiny eccentric- 
ity therefore requires explanation. Since eccentricity 
likely decays exponentially, the current low e, even for 
an initial eccentricity as small as e ~ 0.1, requires a 
circularization time tdrc <J ''"sd/4 — 2.5 Myr. 

It is not straightforward to estimate whether such 
a short tcirc is possible. This is because the compan- 
ion is probably completely convective, with convec- 
tive turnover time scales far longer than the 6.2-h or- 
bital period (see below). This makes the energy trans- 
fer less efficient, and thus circularization time scales 
longer, but it is not clear to w hat e xtent. We will first 
use the prescription of Zahn ( 1989 ), in which the cfh- 
ciency is assumed to decrease linearly with the ratio 
of the convective to orbital timescale, and then dis- 



cuss the prescription of Goldreich & Nicholson (1977) 
and Goodman & Oh ( |1997D , in which the efficiency 
decreases almost quadratically with the timescale ra- 
tio. We should stress that at present it is not clear 
that either prescription is reliable; see Goodman & 



Oh for a discussion. 

Following the formahsm of Zahn ( 1989| ; Eq. [21]), 
we write 

= = 21^1(1 + 1) (^)\ (1) 

tcirc e at tf \ a J 

where ^circ is the circularization timescale, e the eccen- 
tricity, Acirc a dimensionless average of the turbulent 
viscosity weighted by the square of the tidal shear, 
t{ — {MqR^/ LcY/'^ the convective friction time, 
q = Mpsr/Mc the mass ratio, and a the orbital sepa- 
ration. Here, all uncertainty is hidden in the parame- 
ter Acirc- In the prescription of Zahn (1989), it can be 
approximated by Acirc ^ Qm9a^/'^(l + 7^/2,20)-^/'^, 
where a is the mixing length parameter and "q = 
2tf I Porb a measure of the timescale mismatch. 

For a main-sequence star companion with param- 
eters in the ranges listed above, we find t{ ~ 0.5 yr 
and r\ ~ 1500. Assuming AfpgR = 1.4 M q (i.e ., 
a ~ 2i?o) and taking a = 2 (as in Zahn | 1989|| ; 
consistent with the observational constraint derived 



by Verbunt & Phinney |1995| for small 77), one in- 
fers Acirc — 6 X 10^**, and a circularization time in 
the range 5 < icirc <^ 17 Myr, too long to under- 
stand the current small eccentricity. If the star is 
bloated, circularization is much faster, because of the 
very strong dependence on the ratio Rc/a: we find 
<circ — 0.07 Myr for a maximally bloated companion. 

The situation is different for the prescription of 
Goodman & Oh (1997), in which Acirc oc 77"^. Ex- 
trapolating in their Fig. 2, we infer Acirc — 4 x 10~^ 
(note their slightly different definition of 77). Thus, 
the inferred circularization times are two orders of 
magnitude longer than those estimated with the pre- 
scription of Zahn ( 1989 ); if correct, it may be difficult 
to understand how the orbit could have been circu- 
larized even if the companion was maximally bloated. 

As mentioned in §^ the star could have become 
bloated due to th e energy dissip ated by the circu- 
larization proper ( Verbunt 1994 ). It is difficult to 
estimate, however, by how much, as it is not clear 
where, how, and on what time scale the tidal energy 
is dissipated, especially in such a low-mass star. For 
high eccentricity, the tides excite low-order pulsation 
modes in the star, which will be dissipated, either by 
direct viscous damping ( Kochanek 1992 ), or, perhaps 
more likely, by non-linear coupling to higher-degree 
modes and damping of these (Kumar & Goodman 



See Phinney ( L99i ) for a pedestrian derivation. 
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1996). For either case, it appears that for a low-mass, 
(almost) completely convective star (^O.SM©), most 
of the energy will be dumped in the outermost layers. 
These will be heated and may expand, which would 
lead to stronger tidal coupling. 

If the expansion involved the whole star, or a sub- 
stantial fraction of it, most likely the star would still 
be bloated, as the thermal time would be much longer 
than the pulsar characteristic age. Also if only the 
outer layers of the star expanded, however, it seems 
likely the star would still be bloated, as otherwise 
it would be difficult to reduce the eccentricity suf- 
ficiently. This is because the tidal luminosity will 
decrease rapidly with decreasing eccentricity and in- 
creasing periastron distance; if the expanded layers 
shrunk too quickly in response, the circularization 
time would have become long again while the ec- 
centricity was still substantial. One way to verify 
whether the companion is still bloated is to measure 
its temperature. 

5.3. Caveat: Fainter Companion 

There is a probability of about one in five of a 
chance coincidence between PSR B1718— 19 and the 
object in our HST images (§^). If so, the companion 
has to be substantially fainter, with r7iF702W > 27.2 
(§ H). One can go through a similar exercise as in § 4.3 
to constrain the companion properties for this case. 
The result is that the companion needs to have mass 
<^ 0.15 M0, close to the minimum allowed by pulse 
timing. In this case, the circularization time scale 
problem is exacerbated. 

Another possibility is that PSR B 1718- 19 is not 
associated with NGC 6342. Indeed, the dispersion 
measure is only 71 cni^"^ pc, while 130cm~'^pc is ex- 



pected for NGC 6342 ( [Taylor fc Cordes 19931 ). Taken 
at face value, a distance of only 3 kpc is implied. If 
PSR B1718— 19 were in the foreground, again the 
companion would be less luminous and thus less mas- 
sive, even more problematic for circularization. 

6. Conclusions 

In summary, using HST observations, we have de- 
tected a faint star at the position of the unusual eclips- 
ing binary radio pulsar PSR B1718— 19. This faint 
star most likely is the pulsar's companion. We have 
shown that it is difficult to explain the highly circular 
present-day orbit if the companion is a Mq — 0.3 Mq, 
Rc — 0.3 i?Q main-sequence star (unless circulariza- 



tion was not by tidal interaction). If it is a bloated, 
Rc ^ 0.4 i?Q star, circularization may be sufficiently 
rapid, depending on the extent to which the efficiency 
of tidal dissipation is suppressed by the orbital period 
being far shorter than the convective turnover time 
in the star. The system thus provides an interesting 
test-case for tidal-interaction theory. 

A measurement of the color of the companion 
would be the best way to determine whether it is 
bloated or not. For a main-sequence star with Toff = 
3640 K and a bloated star with T^s = 3300 K, one ex- 
pects {R-I,I-J,J-H,H-K)o = (1.0,1.1,0.6,0.2) 
and (1.3, 1.4, 0.6, 0.2), respectively. Radio and optical 
proper-motion studies could settle association of our 
candidate counterpart with PSR B1718— 19, and of 
the system with NGC 6342. 
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